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Data on the hadronic final state of deep-inelastic events at the e±p collider HERA 
are reviewed. Fragmentation properties in the current region extracted from 
charged particle spectra are compared to e + e — and fixed target experiments at 
lower center of mass energies. A measurement of the strong coupling constant 
from integrated jet rates is presented and prospects to extract the gluon density 
in the proton by tagging vector mesons are discussed. Data on the inclusive mean 
transverse energy together with the transverse momentum distribution of single 
charged particles indicate an increased gluon activity in the central rapidity region 
when the phase space is enlarged. The measurement of the forward jet cross- 
section is discussed. Keeping the phase space fixed, the mean transverse energy 
and the mean charged multiplicity for certain rapidity ranges is found to behave 
analogously to pp collisions. 


Introduction 

The large center of mass energy (yfs ss 300 GeV) available at the e ± p collider 
HERA offers a large phase space for production of hadrons in deep-inelastic 
scattering (DIS). This allows the detailed characteristics of the hadronic final 
state to be probed in a previously unexplored kinematic domain. 

In the quark parton model (QPM) the photon scatters elastically on a 
quark freely moving in the proton. This is a good approximation when small 
distances are probed at large momentum transfers Q 2 . The kinematics of the 
scattering process is described by the longitudinal momentum fraction of the 
quark with respect to the proton momentum ( x ) and by the resolution of the 
photon probing the proton structure ( Q 2 ). These variables are defined in a 
general way in Fig.|l] as Lorentz invariants. At HERA, x values down to a few 
10 -5 can be reached and Q 2 ranges from low values where the photon is quasi 
on-shell to very large ones resolving - with the present statistics - distances] 
down to « 10 -16 cm. 

Perturbative QCD plays a key role in the understanding of the complex 
mechanism governing the multi-particle final state. For very short times (cor¬ 
responding to short distances) before and after the hard interaction of a virtual 
photon with a parton in the proton, the partons are far away from their mass 
shell and can be treated as quasi free particles. This first primary scatter- 

“Invited talk given at XVI Int. Conf. on Physics in Collision, Mexico City, June 1996. 
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ing process can be perturbatively described by QCD. Then long range forces 
become more and more important and quarks with very low virtualities form 
bound states leading to the observable hadrons. 

The multi-purpose HERA detectors allow, in addition to the scattered 
electron, a large part of the hadronic final state to be measured. Provided that 
the transition from partons to hadrons can be properly modeled, this opens the 
possibility to study the underlying parton dynamics using exclusive observables 
built on particle spectra, jets or energy depositions. Variables based on the 
hadronic final state can be defined for specific rapidity ranges such that e.g. 
properties of hadrons close to the proton remnant can be investigated. It is 
also possible to isolate the current region where the hard subprocess induced 
by the photon virtuality has taken place and to compare with results obtained 
in e + e _ annihilation. An advantage of HERA is that most measurements can 
be performed at a variable scale. The mean multiplicity of charged particles 
can e.g. be studied as function of W and the strong coupling constant a s can 
be extracted from integrated jet rates in different bins of Q 2 . The phenomenon 
of diffraction in DIS is covered by M. Derrick in these proceedings. 

Particle spectra in the current region 

A natural reference frame to study properties of the hadronic final state is the 
hadronic center of mass system (cms) where the virtual photon and the proton 
have momenta of equal magnitude and point in opposite directions. Within 
the QPM, the photon and the incoming quark collide in this frame head-on and 
no transverse momentum (Bp) is transfered to the struck quark. The Bp of 
the measured hadrons is only produced during the fragmentation of the quark 
and therefore maximally limited to the typical inverse size of a hadron. Higher 
values of Bp can only originate from gluon emissions before and after the hard 
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Figure 2: Mean squared transverse momentum as function of the scaled longitudinal mo¬ 
mentum for charged particles in the cms. HERA data are indicated by the open squares 
(HI) and the closed circles (ZEUS). The data from a p^p fixed target experiment at 
y/~s = 23 GeV (EMC) are given as open circles. The prediction from ARIADNE as a 
QCD model (solid) and the naive quark parton model (dashed) are superimposed as lines. 


interaction. Neglecting gluon radiation the longitudinal momentum (P*) of 
the struck quark is given by W/ 2. Correspondingly the target remnant moves 
in opposite direction with —W/ 2.A natural variable to study the longitudinal 
momenta of the emerging hadrons is therefore Xf = P* / P*m&x = 2P*/W. 

The importance of processes involving QCD radiation is illustrated in 
Fig. I where the (F^ 2 ) versus the xf of charged particles is shown. With¬ 
out allowing for parton radiation the (P^ 2 ) is limited 0.3 GeV 2 (dotted 
line). Much larger (F^? 2 ) are found in the HERA data mil and they are well 
described by a QCD model (solid line). In the target hemisphere (xf < 0), 
where no HERA measurements are available yet, the (Fjf 2 ) is less prominent 
and only little P£ is produced by gluon radiation. This can be understood in 
terms of the extended colour charge distribution supressing small wavelength 
gluon radiation. At lower center of mass energies where less phase space for 
gluon radiation is available, the (Pf 2 ) is even more limited. This can be seen 
e.g. in the data of the EMC collaborations, a / i±p fixed target experiment 
at y/s = 23 GeV. In this kinematic region QCD models do not significantly 
add transverse momenta to the naive expectation from the QPM and have 
moreover difficulties to describe the dataQ. 

Detailed comparisions between e + e _ and e±p collisions are most suitably 
carried out in the Breit frame. In this frame the only non-zero space-like 
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Figure 3: (a) Sketch of a simple photon-quark collision in the Breit frame (b) Phase space 
for the longitudinal and transverse momenta of single particles in this frame. 


component of the photon four-momentum is given by P z = —2xEp with Ep 
the proton beam energy. Within the QPM, the struck quark enters with P z = 
Q/2 and collides on the virtual photon as on a ’brick wall’, from which it 
simply rebounds with P z = —Q/2 (see Fig. ga). The particles scattered in 
the hemisphere of the struck quark (P z < 0) behave like those produced in 
one hemisphere of an e + e~ annihilation event at s =Q 2 . As in e + e _ in this 
hemisphere properties of longitudinal particle momenta should only depend on 
Q. The opposite hemisphere {P z > 0) contains the proton remnant and moves 
with P z = (1 — x)E P = (1 — x)Q/2x (see Fig. g::>). 

The linear rise of (twice) the mean charged particle multiplicity ((n)) with 
log Q for DIS measured at HERi®0 is shown in Fig. ga. For Q > 7 GeV, 
DIS and e + e~ data3 reasonably agree. However, when going to lower Q the 
differences between e + e~ and e±p collisions become more and more distinct. 
A higher multiplicity is observed in e + e _ data taken by MARK1. A clear 
dependence of (n) on x is observed when comparing the HERA data at low 
x with the vp data of a fixed target experiment at FNAlE at x close to 0.1. 
While agreeing at high x with the e + e _ data, in DIS the current hemisphere 
is more and more depopulated at low x where higher order parton emission 
become increasingly important. This is due to an asymmetric migration from 
the current to the target hemisphere, not present in e + e _ annihilation. 

The distribution of the scaled longitudinal momenta x p = P z /P zmax = 
2P z /Q of single charged particles allows fragmentation to be studied in fur¬ 
ther detail. The fragmentation function l/cr to t dcr/dx p , with l/crtot the total 
cross-section in a particular (x,Q 2 ) bin, characterizes the complete process of 
hadron formation including QCD parton evolution and the non-perturbative 
transition from partons to hadrons. It is not calculable in perturbative QCD, 
but can, as the structure function F 2 , be evolved in Q 2 given a measured start¬ 
ing distribution, l/crtot do/dx p rises with decreasing x p and turns over at low 
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Figure 4: Mean multiplicitiy of charged particles (a) and the peak position of the lnl/x p 
distribution (b) as function of Q in e~*~e — , e±p and up collisions. The e±p data are measured 
in the Breit frame and the multiplicity is multiplied by a factor of two. 


x p . This region can be expanded by transforming x p to In l/x p (see Fig.|[i and 
b). At small x p , 1/otot da/dx p is affected by coherence of soft gluons while at 
large x p it can be used to determine a s . 

At small x p the phase space for soft gluon emission is reduced to an angular 
ordered region caused by destructive interference. This suppression of soft 
gluons is predicted by QCD, since the power of gluons to_Lfisolve individual 
colour charges should be diluted with increasing wavelengtlrJ. Assuming that 
the basic features of a distribution for partons are preserved for the measured 
hadrons, this results then in a slower rise of the particle multiplicity and of 
the peak position of the In l/x p distribution with log Q than simply expected 
from the increased longitudinal phase space. Fig. |]b shows the evolution of the 
hil/x p peak position measured for particles having x p > 0. The same linear 
rise with logQ is observed in e + e _ datalu which provides a hint of a universal 
behaviour of quark fragmentation. As the data become more precise, expected 
differences between e“*"e _ and DIS fragmentation should become manifest due 
to e.g. different quark flavour composition. The modified leading logarithm 
approximation (MLLA) of QCD is able to describe this evolution down to very 
low Q (solid line). The coherent evolution is clearly favoured by the data. It 
can, however, not be excluded that a similar attenuation of the slope could 
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Figure 5: Distribution of the scaled longitudinal momentum x p = 2P Z IQ for charged par¬ 
ticles in the Breit frame in two (#,Q 2 ) bins measured by ZEUS (b). The same data after 
transforming x p to \nl/x p (a). 


also be caused by effects during the hadronisation 


phascEl. 


Xlic evolution of 1/utot da/dx p in different bins of x p with Q 2 is shown in 
Fig. J. With increasing Q 2 , 1/crtot da/dx p increases at low x p and decreases 
at high x p . When more phase space is available more gluons can be emitted 
which results in a softening of the particle spectra. Therefore more particles at 
low x p and less at large x p are observed. These scaling violations can be treated 
in an analogous way as in the proton structure function analysi£3 to extract 
the strong coupling constant a s . Such a measurement has been performed liJ 
using LEP data together with lower energy e + e _ data. At HERA this method 
could be applied in one single experiment. However, to avoid the region were 
migrations between the current and target hemisphere become important, one 
has to restrict the analysis to Q 2 > 100 GeV 2 . Scaling violations are observed, 
but more data is needed to get the necessary precision at large x and Q 2 for a 
competitive extraction of a s . Since the scaling violation in the fragmentation 
function are only of logarithmic nature, the achievable error on a s will, even 
with high statistics, not be competitive with the present error on the world 


averagaU. 
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Figure 6: 

l/fltot da/dxpin the Breit 
frame for different bins of 
the scaled longitudinal mo¬ 
mentum x. p as function of 
Q 2 . The data have been 
measured by ZEUS. 


Determination of the strong coupling constant using jets 

The high values of W reachable at HERA allow clean jets to be observed in 
the hadronic final state. Events with two jets in addition to the one close to 
the beam carrying the remnants of the proton (2+1 jet events) are produced 
in leading order of a s by processes either induced by a gluon in the proton 
(Fig. 0a) or by a quark radiating a gluon before or after interacting with the 
photon (Fig. 00- The cross-section of both processes is proportional to a s . 
For the gluon initiated processes in addition precise knowledge of the gluon 
density in the proton g(x,Q 2 ) is required to calculate the cross-section. The 
sensitivity to both a s and g(x, Q 2 ) can be used to simultaneously extract 
these quantities. Such an analysis is however theoretically and experimentally 
difficult. The analysis is therefore split in two kinematic regions. At low 
x and Q 2 the gluon initiated processes dominate. This can be used to extract 
g(x,Q 2 ) by assuming a value of a s and statistically subtracting the quark 
induced processes. At large x and Q 2 the parton densities are well known and 
the proportionality coefficients between_tlre 2+1 jet rate and a s are given by 
next-to-leading order QCD calculations!^. The measurement of the ratio of 
2+1 jet events to all events at some value of Q allows therefore a s (Q) to be 
extracted. By evaluating a s (Q) in different Q bins the running of the strong 
coupling constant can be tested in one experiment. 

For the a s measurement only the kinematic range of high Q is considered 
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Figure 7: Feynman diagrams for the production of 2 + 1 jet events to first order of a s in 
e^p-collisions. 

to reduce the uncertainties from hadronisation. To avoid the region where 
higher order gluon emission become important, only events at high x are 
considered, and also, the experimentally difficult, very forward region is ex¬ 
cluded. Jets are_experimentally and theoretically defined by Hie JADE clus^ 
tering algorithnlilThree values of a s (Q) extracted by the HI rill and ZEUScil 
collaboration are shown in Fig. ||. For increasing values of Q, a s (Q) decreases 
as predicted by the renormalisation group equation shown for A^jg between 
100 and 300 MeV (solid lines). An extrapolation to a s (M z ) yields: 

HI 93 : a s {M z ) = 0.123 ± 0.012(stat) ± 0.013(syst.) 

ZEUS 94 : a s {M z ) = 0.117 ± 0.005(stat) + 0.005(exp.) ± 0.007(th.) 

which is consistent with othemvalues obtained from a large variety of different 
processes as shown Fig. || (se£3 for references). The agreement found in e + e _ , 
e±p and pp collision for reactions at very different scales is an important and 
successful test of QCD. The error on the HERA measurement are even with 
the present limited statistics already competitive. 

The dominant experimental errors stem from the uncertainty on the energy 
scale for hadrons, the model dependence of assigning the jets to partons and the 
phenomenological description of the hadronisation process. Furthermore the 
renormalisation scale and the choice of the input parton density systematically 
influences the result. For consistency it would be necessary to simultaneously 
fit a s used in the parton densities, but this effect is negligible in the phase space 
region where the HERA analysis has been performedEj. Recently a new next- 
to-leading order calculatioio became available allowing to analyze the data 
using all available jet definition schemes and to impose any acceptance cut on 
the final state particles. Moreover different renormalisation and factorisation 
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Figure 8: Left: Values and total error of a s (Mz) from various processes. The solid line 
indicates the world average and the band its total error. Right: a s (Q) from HERA (open 
symbols) and other processes with increasing Q (closed circles): Tx and ohad/^tot, event 
shapes and r^adron/riepton i n e+e - • The predictions of the renormalisation group equation 
for three values of A^^- are superimposed as lines. 

scales (e.g. Pf instead of Q) can be tested. These calculations will allow 
the theoretical systematic errors to be better assessed. One of the insights 
gained from these studies is that the como or fcj l 3 jet algorithm seem to 
be better suited for precision QCD tests than the JADE scheme, since the 
variation of the .cross section when going from leading to next-to-leading order 
is much smalle£3. Moreover, these jet algorithms are less dependent on the 
choice of the hard scale and allow a larger phase space to be covered. The 
ZEUS collaboration has reanalyzed their datacS using the kr algorithm. The 
preliminary values of a s (Q) obtained in three bins of Q are showrf] in Fig. || 
They are consistent with the results obtained with the JADE algorithm. 


On the determination of the gluon density in the proton 

The dynamics of the strong interaction at high momentum scales, correspond¬ 
ing to short distances, is successfully described by perturbative QCD. A vital 
ingredient of any QCD prediction for a reaction where a nucleon is involved is 
an assumption about the parton distribution in the nucleon needed to calculate 
an observable in terms of basic QCD subprocesses at the partonic level. This 
parton density function f{x,Q 2 ) is the probability to find a parton within a 

^Only the statistical errors are included in the figure. 
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Figure 9: Leading Feynman diagram for inelastic J/tp production (a) and a QCD inspired 
model for elastic J/tp production (b) in 'ip-collisions. 


nucleon carrying a fraction x of its momentum when probed by a particle with 
virtuality Q 2 . Once f(x,Q 2 ) has been determined in one reaction it can be 
used for any other prediction. 

In deep-inelastic scattering the collinear singularities encountered when 
calculating the hard subprocess to order a s are absorbed in the parton den¬ 
sities. The parton densities are usually evolved by the DGLAP equation^, 
presently known in leading and next-to-leading order, which describe the split¬ 
ting of a parton. Once the parton densities have been measured at some 
Q 2 value they can be evolved to any point in the kinematic plane. The 
parametrisation of the input distribution is obtained by a global fit to a large 
variety of data. Two different strategies haye been adopted. Martin, Roberts 
and Stirlingp and the CTEQ collaboration!!] parametrise the starting distribu¬ 
tion at some Q q (1 — 5 GeV 2 ) and fit then the parameters to data by -evolving 
/( x, Qq). A dynamical approach is pursued by Gliick, Reya and VogtJ where 
only valence like parton densities with f(x = 0, Q q) = 0 are assumed at a very 
low input scale Qq = 0.34 GeV 2 . The gluon and the sea quark distributions 
are then generated radiatively. This procedure predicted the steep rise of the 
structure function [] with decreasing x and moreover the slowing down of 
this rise towards lower Q 2 . 

The parton distributions are well constrained for large x where a lot of data 
are available. However, in the small x regime where gluons play a crucial role, 
they are not well known. HERA has therefore the unique possibility to investi¬ 
gate the distribution of the gluon density in the proton. The scaling violation 

c see p.14-15 for more details and references. 
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Figure 10: Cross-section for inelastic (a) and elastic (b) J/tp production as function of W lp for 
HERA and fixed target photoproduction data. Predictions of the Ryskin model for different 
parameterisations of the gluon density are superimposed. The solid line below the data 
indicates a simple parametrisation a-y P ~ l+) p for S = 0.32. 


of the proton structure function F 2 offers an indirect way to get a handle on the 
gluon density. This measurement was made possible by the extremely success¬ 
ful description of F 2 by the DGLAP equations and has been presented by M. 
Kasemann at this conference. An alternative way is to directly tag the quark 
box formed in a photon-gluon collision either by identifying heavy quarks, e.g. 
charm via D* or J/ip mesons, or by measuring inclusively the 2+1 jet cross 
section at small x. This method has, compared to heavy flavour production, 
a much larger cross-section, but has the disadvantage that the background 
from quark initiated processes has to be subtracted statistically. Moreover, 
the reach to small values of x g / p =x{l + s/Q 2 ) (see Fig.|?] for notation) is re¬ 
stricted, since large invariant jets masses (s > 100 GeV 2 ) are experimentally 
required to define clean jets. An extraction of the gluon density in leading 
order using the inclusive 2+1 jet rate achieved by the HI collaboiatioitj has 
already been presented on the last ’Physics in Collision’ conferencfa 

Gluon initiated processes can be unambiguously pinned down by selecting 
events where a J/ip is produced. The leptonic decay modes (J/ip —► e + e _ and 
J/ijj —> pi + pi~) of this meson provide a clean experimental signature. In the 
elastic channels only two leptons are found in the detector, in the inelastic 
channel there are additional particles (see Fig. []). Since the charm quark mass 
is large enough to provide a hard scale, this measurement can be performed at 
very low Q 2 where the cross-section is large. Besides the requirement of at least 
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Figure 11: Total and vector meson photoproduction cross-section as function of W 1P mea¬ 
sured at HERA and by fixed target experiments. Superimposed is a simple parametrisation 
cT-yp ~ Wyp for S = 0.32 (light vector mesons) and 8 = 0.6 (J/ip meson). 


one additional particle in the final state, the inelastic channel can be inclusively 
defined using the inelasticity variable z = P ■ Pj/^/{P ■ q) < 0.8 where P ( q ) is 
the incoming proton (photon) momentum and Pj/tf, the four-momentum of the 
J/'i/’. To avoid the phase space region where multiple emissions of gluons be¬ 
come more and more important and where a resummation would be required, 
the P£ of the J/^> shoulcLexceed 2 GeV. In this region a QCD calculation^ can 
predict the cross-sectioiEn of J/^> production in next-to-leading order. The de¬ 
pendence of the cross-section on the photon-proton invariant mass ( W lv ) and 
the overall normalisation is correctly described (see Fig. |ic|a). In this calcula¬ 
tion the J/ip was assumed to be produced in a colour singlet state, no room is 
left for an additional component of a colour octet as required by pp dat£3 at 
high energy. However, excluding the region of small P£, significantly reduces 
the sensitivity to the small-a; behaviour of the gluon density as can be seen 
in Fig. 10 1 where cross-section calculations using various parametrisations are 
compared to data. 

Previous measurements^ of elastically produced J/'ip mesons revealed an 
exponentially falling distribution of the momentum transfer between the in- 
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coining and outgoing proton as expected from diffractive processes and a trans¬ 
verse polarisation in agreement with s-channel helicity conservation. However, 
the W rp dependence of the elastic cross-section a e i measured over a wide range 
by fixed target experimentO and at HERAc9 exhibits a much steeper rise (see 
Fig. fiojb) than predicted by models based on diffractive mechanisms like e.g. 
the one by Donnachie and Landshofcd which assumes a soft interaction in 
a Regge form where no quantum numbers are exchanged. This model relies 
on the observation that any total hadronic cross-section follows a power law 
with the energy available in the center of mass s like cr ~ s A with A ~ 0.08. 
If one parameterizes a e i in the region 30 < W lp < 260 GeV as a e i ~ W pp 
one finds S = 0.6 ± 0.12 rather than S = 0.32 as predicted by the Regge- 
inspired modeL-This significant rise can be better described by Ryskin’s QCD 
inspired modeea where the J/i/> couples via two gluons to the proton. To 
leading order it predicts that a e i ~ [a 3 {Q 2 ) xg(x,Q 2 )] 2 where Q 2 = 1/4 M 2 ^ 
and x = M 2 /^/W 2 p . The quadratic dependence on g(x, Q 2 ), makes the elastic 
channel of J/ij) photoproduction a sensitive discriminator for the slope of the 
gluon distribution a e i ~ [xgix,Q 2 )Y ~ x~ 2X ~ W^ x /Q 2 . This model, includ¬ 
ing higher ordeivcorrectionsiil agrees with the data when using the MRSA’ 
parametrisatioiEJ with A = 0.17. The parametrisation of GRV tends to give a 
steeper W 1P dependence. Since phenomenological assumptions concerning the 
two gluon interaction have to be made, the W lp dependence of a e i, rather than 
the absolute normalization, is a more reliable probe of g(x,Q 2 ). When com¬ 
paring to HERA data, the model prediction is therefore normalized to fixed 
target data. A direct extraction of g(x, Q 2 ) seems feasible once experimental 
and theoretical problems are better understood. 


The steep rise of a e i with increasing W lp can be contrasted with the mod¬ 
erate increase of the total photoproduction cross-section and the cross-section 
of the light vector mesons p, uj and ip which are in excellent agreement with 
model predictions where only soft interactions are assumed (see Fig. 0>- The 
failure of this picture to describe elastic J/ip production together with the suc¬ 
cess of Ryskin’s model gives additional confidence in the potential of QCD to 
describe processes where the scales involved are still surprisingly small. The 
production of light and heavy vector mesons at HERA will allow to gain a 
better understanding of the transition from a region governed by soft interac¬ 
tions to a domain where perturbative QCD turns on. First indications for a 
similar transition when studying light-vector meson production with increasing 
Q 2 have also been observed at HERAe3. 
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Figure 12: Diagram of an e±p collision at low x. 


Deep-inelastic scattering at small x 

At low x the simple picture of deep-inelastic scattering as a process where a 
virtual photon interacts instantaneously with a point-like parton freely moving 
in the proton has to be modified. Since the probability that a parton radiates 
becomes increasingly high at low x, the parton struck by the photon origi¬ 
nates most likely from a cascade initiated by a parton with high longitudinal 
momentum. The phase space for gluon radiation is enlarged with decreasing 
x. Such an interaction is illustrated in Fig.[L2]. The gluon is the driving force 
behind the parton cascade. The transitions g —> gg and g —> qq at each point 
in the ladder can be approximated by the DGLAP equations. By resumming 
the LO and NLO (a s InQ 2 ) terms they predict the Q 2 evolution of a parton 
known to be point-like at some given scale Q o to the region where the inter¬ 
action with the photon takes place. They describe the change in the parton 
density with varying spatial resolution (A ~ 1/Q) of the probe. One of the 
assumption^ necessary to derive the DGLAP equations is a strongly ordered 
configuration in the parton virtualities along the ladder connecting the soft 
proton constituents to the hard subprocess. This leads to a suppression of the 
available phase space for gluon radiation. 

When the (a s \nl/x) n terms become large, they have to be taken into ac¬ 
count, e.g. by the resummation accomplished by the BFKL equations!!. In a 
physical gauge, these terms correspond to an n-rung ladder diagram in which 
gluon emissions are ordered in longitudinal momentum. The strong ordering 

'Tor a digestible introduction to low-.r physics see A. MartirH. 
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Figure 13: Mean transverse energy in —0.5 < r/* < 0.5 measured by HI and ZEUS (closed 
symbols) in the cms for Q 2 = 14 GeV 2 as function of x. Superimposed are QCD model 
predictions (lines) for hadrons (a) and partons (b) by QCD models and an analytic BFKL 
calculation (open circles) carried out on parton level only. 


of the transverse momenta is replaced by a diffusion pattern as one proceeds 
along the gluon chain. The BFKL equations describe how a particular high 
momentum parton in the proton is dressed by a cloud of gluons at low x lo¬ 
calized in a fixed transverse spatial region of the proton. At very low x it is 
expected that many gluons coexist and will no longer act as free partons, but 
interact with each other. This ’saturation’ regime is characterized by an equi¬ 
librium of gluon emission and absorption. A new QCD regime is reached where 
individual parton-parton interactions are weak, but where the field strength 
becomes - due to the number of partons - so strong that perturbation theory 
is not reliable. It is not clear if in the HERA regime such effects could be 
observed. 

Experimentally, it is interesting to measure observables sensitive to the 
underlying parton dynamics.Studying the probe after the interaction is the 
classical approach to get information on the structure of the target. The mea¬ 
surement of the proton structure function F -2 has the benefit that it can be 
directly compared to analytical QCD calculations. The rise of Fi with de¬ 
creasing x observed at HERA3 appears to be well described by the DGLAP 
equations. Even at the lowest values of x the BFKL terms are not needed to 
obtain a satisfactory description. 
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Figure 14: Transverse momentum distribution of single charged particles in 1.5 < rj* < 2.5 
measured by HI in the cms for fixed Q 2 = 14 GeV 2 and two x bins. Superimposed as lines 
are QCD model predictions. 

However, F -2 might be too inclusive to reveal the potentially small resum¬ 
mation effects. A better indicator might be provided by observables based 
on the hadronic final state emerging from the parton cascade, since they can 
be chosen such that the differences between the two complementary approx¬ 
imations are enhanced. However, predictions have then mainly to be based 
on Monte Carlo event generators simulating the detailed properties of the 
hadronic final state by including the hard subprocess, QCD radiation effects 
and the transition of partons to hadrons. 

The LEPTO Monte Carl(@ incorporates the QCD matrix element to first 
order in a s and approximately accounts for higher order parton emissions 
by the concept of parton showerS based on the DGLAP equations (MEPS 
model). The non-perturbative transition of the partonic. final state to the 
observahlc hadrons is modeled by the LUND string modecEl as incorporated in 
JETSETE3. ARIADNfa differs only in the treatment of the hard subprocess 
and the QCD^ shower evolution. Gluon emissions are treated by the colour 
dipole modcLj assuming a chain of independent radiating dipoles spanned by 
colour connected partons. TheJirst emission in the cascade is corrected by 
the matrix element to first ordefiH A phenomenological suppression of gluons 
with short wavelengths (having high transverse momenta) only resolving part 
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of the extended colour charge distribution of the proton remnant has been 
introduced for DIS. This leads to a more pronounced population of partons 
towards the proton remnant than observed in LEPTO where the phase space 
is controlled by the strong ordering of the gluon virtualities and the ratio of 
the parton density function at each branching point. It has been argued that 
the partonic state as encountered i n AR TADNF, is more closely related to the 
one expected from a BFKL scenaridEOH. Both Monte Carlo models are able to 
give a consistent picturmofjthe event topology and provide a fair description 
of basic event propertie££(Ej. 

Using the hadronic final state, the presence of parton activity can be ex¬ 
perimentally probed from the region close to the photon quark interaction 
down to (pseudo-)rapidities defined as ry* = — In tan ( 0 / 2 ) where 0 is the polar 
angle between a particle and the proton direction^. In the cms ?y* sa — 1 can be 
reached beyond which the acceptance of the two main HERA detectors ends 
(see also Fig©. An inclusive measure is the mean transverse energy (( E/)) 
in the central region —0.5 < ry* < 0.5. Due to the increasing phase space with 
decreasing x, gluons are more a bu ndant and (E/) is expected to rise. This 
behaviour is seen in the dataOHS as is illustrated in Figjl3^i where the (El/) 
is shown for fixed Q 2 = 14 GeV 2 . Both QCD models are able to describe 
the data. LEPTO, however, has to produce 60 — 80% of the (E/) during the 
hadronization phase to compensate for the different x dependence of the (E/) 
seen at the parton level as is illustrated in Figjl^. This could only be achieved 
by introducing two new non-perturbative production mechanisms for (E/). In 
events where a sea-quark is involved in the hard subprocess, its partner is not 
- as in previous versions - simply rearranged to a meson or a baryon within the 
proton remnant, but is used to stretch a string to one of the quarks remaining 
untouched in the proton (SOT). The second new ingredient is the assumption 
of a soft colour interactioiEll(SCI) which changes only the colour configuration 
of the partonic system while leaving the colour field of the proton. This colour 
rotation offers, by producing a colour singlet in the final state, a possible ex¬ 
planation of the rapidity gap events observed at HERA and at the same time 
(due to longer strings caused by the colour rearrangement) the transverse en¬ 
ergy needed to match the data. Without these two new features a constant 
amount of (E/) is generated during the hadronisation and the data are not de¬ 
scribed (see Fig.p~3|). Although ARIADNE uses the same hadronisation model 
as LEPTO, a constant amount of only 30-40%xtf the (E/) is produced during 
hadronisation. An analytical BFKL calculation^, only carried out for partons, 
predicts the same x dependence as seen in the data, but lies significantly above 
the corresponding ARIADNE curve. 

6 The direction of the proton points along the negative z-axis in the cms. 
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Figure 15: Forward jet cross-section as function of x measured by HI. Superimposed are 
QCD model predictions for hadrons (a) and partons (b). An analytic BFKL calculation 
(open squares) and a NLO QCD calculation (open circles) are superimposed. 


The hard tail of the distribution of single charged particles has been 
suggested^ as a better tool to disentangle non-perturbative hadronisation ef¬ 
fects from parton activity as described by perturbative QCD. This observable 
is more directly related to hard parton emissions and cannot be mimicked by a 
cumulative effect from many soft particle emissions during hadronization. The 
measured steeply falling distribution in the interval 1.5 < rj* < 2.5 gets 
harder for decreasing x (see Fig.^fj). While at large x LEPTO and ARIADNE 
are both able to describe the data, LEPTO significantly falls below the data 
at low x. This indicates that more partons are produced than expected from 
models with suppressed parton radiation based on the DGLAP parton showers. 
Similar, but less conclusive, resultfj have been obtained by using the inclusive 
transverse energy distribution in —0.5 < r/* < 0.5. 

Another signature designed to enhance BFKL effects is the production 
of ’forward jets’ E3 characterized by a squared transverse momentum (k\.) of 
the order of the photon virtuality Q 2 and Xj e t = E- ]et /E p as large and x as 
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Figure 16: Energy flow at fixed W = 180 GeV 2 in the hadronic center of mass for different 
Q 2 bins ranging from 0 to 520 GeV 2 and mean transverse energy for —0.5 < rf < 0.5 and 
2 < Tj* < 3 in function of Q 2 . 


small as kinematically possible. E - ]e t (E p ) denotes the energy of the forward 
jet (proton). The first requirement suppresses the strongly ordered DGLAP 
evolution. Asking for large Xjet/x, the forward jet can be separated by a large 
rapidity interval from the struck quark such that the phase space of parton 
emission between the two is amplified. For this kinematical configuration the 
a s In x- ]e t /x terms are expected to become so large that their resummation 
should lead to a sizable increase of the forward jet cross-section. 


A fast rise of the cross-section for forward jets defined by kx > 3.5 GeW. 
0.5 < kx/Q 2 < 2 and lEjet > 0.035 with decreasing x is seen in the HI dataS 
corrected to hadron level (see Figjisla). ARIADNE falls slightly below the 
data, but shows the same trend towards small x. Hadronisation effects are well 
below 20% such that the same x dependence is also seen on parton level. While 
LEPTO predicts - as expected - a similar moderate increase of the forward jet 
cross-section as a full NLO QCD calculation!!] at the parton level, it is also 
able to describe the data by producing a large fraction of forward jets (up to 
80%) in the hadronisation phase. Without assuming soft colour interaction the 
hadronisation corrections are similar to the ones obtained in ARIADNE, but 
then the chita are not described (dotted line in Fig.|l5|a). An analytical BFKL 
calculation!))] exhibits a much faster rise with 1 /x compared to the ARIADNE 
result on parton level. However, the authors point out that several effects 
which might lower the prediction have not been taken into account. 
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Figure 17: Mean multiplicity of charged particles in a given (pseudo-)rapidity region y(rf) 
and ( E? r ) in a central rapidity region measured in er p and pp collsions as function of W. 
The solid (dotted) line indicates the behaviour of the mean charged particle multiplicity in 
pp (e+e — ) collsions. 

Deep-inelastic scattering as hadron-hadron scattering 

Strictly speaking, the physical picture of DIS as an interaction of a photon 
with a parton freely moving in the proton is only valid if the proton is very 
fast. Only in this reference frame is the interaction time of the virtual photon 
with the parton short compared to the interactions regularly taking place in 
the proton such that the photon instantaneously probes the proton content. In 
the reference frame where the proton is at rest and the photon is fast, it is more 
appropriate to imagine that the photon fluctuates into a quark anti-quark pair 
before interacting with the proton. If the fluctuation time is long compared to 
the interaction time, a parton cascade can develop and a complicated object is 
formed which strongly interacts with matter. Using the uncertainty principle 
the fluctuation tim<££l can be calculated to r 7 = l/(m p x) where m p is the 
proton mass. For a typical x value of 10 -3 , the distance before the photon 
interacts with the target is cr 7 ~ 200 fm which is large compared to the 
proton radius of about 1 fm. 

If one chooses suitable observables, one can therefore find in DIS at low- 
x characteristic properties of hadron-hadron collisions. Changing thmphoton 
virtuality from « 0 to about 520 GeV 2 , the transverse energy flovEJ in the 
hadronic cms for fixed W = 180 GeV varies only in the region where the hard 
subprocess has taken place (see Fig. 0§i). The (Ej) in 2 < rj* < 3 indeed 
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grows with Q 2 as is depicted in Figjl^b. Surprisingly, this rise only starts 
at about 10 GeV 2 . The measured (Ef) produced in collisions with on-shell 
photons, which cannot be treated perturbatively, has the same magnitude as 
observed in a high Q 2 region where first calculations have been very promising. 
In the central rapidity region of —0.5 < rf < 0.5 however, (Elf) is only little 
influenced by the virtuality of the photon. The impact of the photon virtuality 
penetrates only A rf ss 1.5. This is typical of the short range correlation in in¬ 
elastic hadron-hadron collisions. When comparing the (Ef) to hadron-hadron 
collisions, it scales with the center of mass energy of the collision indepen¬ 
dent of the nature of the incoming particles. This is shown in Figjl7[ It is 
interesting to note that also the mean charged particle multiplicity exhibits a 
similar dependence on W. Also here the DIS data from HERAO interpolate 
- despite the slightly different rapidity regions - between previous fixed target 
data and hadron-hadron collisions at much higher energy (see Fig.|l7|). This 
suggests a universal behaviour of the dynamics in the central rapidity region. 
In e + e^ collisions a steeper dependence on W is seen for these rapidity regions 
(dotted line in Fig.pT). 

Since DIS data reveal some properties of soft hadron-hadron interactions 
there is some hope that the understanding of DIS at low x by perturbative 
QCD can give a handle on the understanding of the underlying dynamics of 
strong interactions. If the first attempts to understand in perturbative QCD 
some low-a; HERA data are successful, a tool might become available to provide 
a universal description of strongly interacting particles from an perturbative 
theory of asymptotically free partons at short distances to the theory of the 
confined hadrons at long distances. 

Acknowledgments 

It is a pleasure to thank my colleagues from the HI and ZEUS collaboration for 
providing me with the latest data and information. I have profited from many 
stimulating discussions with M. Kuhlen and E. DeWolf and I gratefully appre¬ 
ciate their careful reading of the manuscript. I deeply acknowledge the support 
of the physics coordinators of the HI and ZEUS collaboration, J. Dainton and 
R. Nania, and would like to thank them for their helpful comments. 

References 

1. HI Collab., S. Aid et al., Phys. Lett. B353 (1995) 578; ZEUS Collab., 
Prelim, results, see M. Kasemann in these proceedings. 

2. HI Collab., I. Abt et al., Z. Phys C 63 (1994) 377. 

3. ZEUS Collab., M.Derrick et al. Z. Phys. C70 (1996) 1. 


21 


4. EMC Collab., A. Arnedodo et al., Z. Phys. C36 (1987) 527. 

5. N. Magnussen et al., Proc. of the HERA workshop 1991, ed. G. Ingel- 
mann and W. Buchmuller, Vol. 3, Hamburg 1992, p. 1167. 

6. HI Collab., S. Aid et al., Nucl. Phys. B445 (1995) 3. 

7. ZEUS Collab., M.Derrick et al. Z. Phys. C67 (1995) 93. 

8. MARK1 Collab., J.L. Siegrist et al., Phys. Rev. D26 (1986) 969; OPAL 
Collab., P.D. Acton et al., Z. Phys C53 (1992) 539; PLUTO Collab., C. 
Berger et al., Phys. Lett. B95 (1980) 313; TASSO Collab., W. Braun¬ 
schweig et al., Z. Phys. C45 (1989) 193. 

9. FNAL Collab., M. Derrick et al., Phys. Lett. B91 (1980) 470. 

10. Y. L. Dokshitzer et al., Basics of perturbative QCD, Editions Frontieres, 
Gif-sur-Yvette, 1991. 

11. TASSO Collab., W. Braunschweig et al., Z. Phys. C67 (1995) 81; Z. 
Phys. C22 (1984) 307; OPAL Collab., M. Akraway et al., Phys.Lett. 
B247 (1990) 617; TOPAZ Collab., R, Itoh et al., Phys. Lett" B345 
(1995) 335. 

12. HI Collab., S. Aid et al., Nucl. Phys. B445 (1995) 3. 

13. R.D. Ball and S. Forte, Phys.Lett. B358 (1995) 365. 

14. DELPHI Collab., P. Aarnio et al., Phys. Lett. B311 (1993) 408. 

15. ZEUS Collab., M. Derrick et al., Z. Phys. C67 (1995) 93. 

16. N. Brook, Proc. of the Int. Workshop on DIS and Related Phenomena, 
ed. G. D’Agostini and A. Nigro, Rome, 1996 and ZEUS Collab., PA02- 
034 contr. to ICHEP’96, Warsaw, Poland, July 1996. 

17. D. Graudenz, Proc. of ’Future Physics at HERA’, ed. G. Ingelmann, A. 
De Roeck and R. Klanner, Hamburg 1996, p. 533. 

18. D. Graudenz, Phys. Lett. B256 (1991) 518; Phys. Rev. D49 (1994) 
3291; T. Brodkorb and J. G. Korner, Z. Phys. C54 (1992) 141. 

19. JADE Collab., W. Bartel et al., Z. Phys. C33 (1986) 23. 

20. HI Collab., T. Ahmed et al., Phys. Lett. B346 (1995) 415. 

21. ZEUS Collab., M. Derrick et al., Phys. Lett. B363 (1995) 201. 

22. Particle Data Group, R.M. Barnett et al., Phys. Rev. D54 (1996). 

23. J. Chyla and J. Rames, Proc. of ’Future Physics at HERA’, ed. G. 
Ingelmann, A. De Roeck and R. Klanner, Hamburg 1996, 529. 

24. E. Mirkes and D. Zeppenfeld, Phys. Lett. B380 (1996) 205. 

25. J. Huth et al., in Proc. Summer Study on HEP, Snowmass, 1990. 

26. S. Catani, Y. Dokshitzer and B. Webber, Phys. Lett. B285 (1992) 291. 

27. E. Mirkes and D. Zeppenfeld, TTP 96-10, MADPH-96-935, invited talk 
given at the Cracow Epiphany Conference on Proton Structure, Jan 1996. 

28. T. Trefzger, Proc. of the Int. Workshop on DIS and Related Phenomena, 
ed. G. D’Agostini and A. Nigro, Rome, 1996. 


22 



30. 


31. 


34. 

35. 

36. 


29. Yu. L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641; V.N. Gribov and 

L. N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438 and 675; G. Altarelli 
and G. Parisi, Nucl. Phys. 126 (1977) 297. 

A. D. Martin, R.G. Roberts and W.J.Stirling, Phys. Rev. D 50 (1994) 
6734 and Phys. Lett. B354 (1995) 155. 

CTEQ Collab., H.L. Lai et al., Phys. Rev. D51 (1995) 4763. 

32. M. Gluck, E. Reya and A. Vogt, Z. Phys. C67 (1995) 433. 

33. HI Collab., S. Aid et al., Nucl. Phys. B449 (1995) 3. 

U. Gorlach, DESY 95-207, Hamburg 1995. 

M. Kramer et al., Phys. Lett. B348 (1995) 657; M. Kramer, Nucl. Phys. 
B459 (199) 3. 

HI Collab., S. Aid et al., PA02-085 contr. to ICHEP’96, Warsaw, July 
1996; ZEUS Collab., M. Derrick et al., Phys. Lett. B350 (1995) 120 and 
PA02-047 contr. to ICHEP’96, Warsaw, July 1996; 

37. M. Cacciari and M. Kramer, Phys. Rev. Lett.76 (1996) 4128. 

38. HI Collab., S. Aid et al., Nucl.Phys. B472 (1996) 3. 

39. E687 Collab., P.L. Frabetti et al., Phys. Lett. B316 (1993) 197; E401 
Collab., M. Binkley et al., Phys. Rev. Lett. 48 (1982) 73; FTPS Collab., 

B. H. Denby et al., Phys. Rev. Lett. 52 (1884) 795. 

40. HI Collab., S. Aid et al., PA02-085 contr. to ICHEP’96, Warsaw, July 
1996. ZEUS Collab., M. Derrick et al., Phys. Lett. B350 (1995) 120 and 
PA02-047 contr. to ICHEP’96, Warsaw, July 1996. 

41. A. Donnachie and P.V. Landshoff, Phys. Lett. B348 (1995) 213; ibid. 
B296 (1992) 227. 

42. M.G. Ryskin, Z. Phys. C57 (1993) 89._ 


43. M.G. Ryskin et al., DTP/95/96 and hep-ph/9511228 


44. A.D. Martin, W. J. Stirling and R. G. Roberts Phys.Rev.D51 (1995) 4756. 
DTP/95/96 and |hcp-ph/951122Sl . 

45. ZEUS Collab., M. Derrick et al., Phys. Lett. B356 (1995) 601. 

46. A. Martin, Lectures given at the XXI Int. Meeting on Fundamental 
Physics, Miraflores de la Sierra, Madrid, 1993 and DTP/93/66, August 
1993. 


47. E.A. Kuraev, L.N. Lipatov and V.S. Fadin, Sov. Phys. JETP 45 (1972) 
199; Y.Y. Balitsky and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 282. 

48. HI Collab., S. Aid et al., Nucl. Phys. B470 (1996) 3; ZEUS Collab., M. 
Derrick et al., Z. Phys. C69 (1996) 607. 

49. G. Ingelman, in Proc. HERA workshop, Eds. W. Biichmiiller and G. 
Ingelman, Hamburg (1991) Vol. 3, 1366; G. Ingelman, A. Edin and J. 
Rathsman, DESY-96057, Hamburg, April 1996. 

50. B. Andersson et al., Phys. Rep. 97 (1983) 31. 


23 





51. A. Edin, G. Ingelman and J. Rathsman, DESY 96-060, Hamburg, April 
1996. 

52. M. Bengtsson and T. Sjostrand, Z. Phys. C37 (1998) 465; M. Bengtsson, 
G. Ingelman and T. Sjostrand, Proc. of the HERA workshop 1987, ed. 
R.D. Peccei, DESY, Hamburg Vol. 1 (1988) 149. 

53. T. Sjostrand, Computer Phys. Comm. 39 (1986) 347; T. Sjostrand and 
M. Bengtsson, ibid. 43 (1987) 367; T. Sjostrand, CERN-TH-6488 (1992) 
and CERN-TH 7112 Dec. 1993. 

54. L. Lonnblad, Computer Phys. Comm. 71 (1992) 15. 

55. G. Gustafson and U. Peterson, Nucl. Phys. B306 (1988); G. Gustafson, 
L. Lonnblad and U. Peterson, Z. Phys. C43 (1989) 625. 

56. L. Lonnblad, Z. Phys. C65 (1995) 285. 

57. A. H. Mueller, Nucl. Phys. B 415 (1994) 373. 

58. J. Rathsman, SLAC-PUB-7344, October 1996. 

59. T. Carli, Proc. of the Int. Workshop on DIS and Related Phenomena, 
ed. G. D’Agostini and A. Nigro, Rome, 1996. 

60. N. Brook et al., Proc. of ’Future Physics at HERA’, ed. G. Ingelmann, 
A. De Roeck and R. Klanner, Hamburg 1996, 613. 

61. HI Collab., S. Aid et ah, Phys. Lett. B356 (1995) 118. 

62. HI Collab. contrib. PA02-073 to the 28th ICHEP’96, Warsaw, July 
1996. 

63. N. Pavel, Proc. of the Int. Workshop on DIS and Related Phenomena, 
ed. G. D’Agostini and A. Nigro, Rome, 1996. 

64. K. Golec-Biernat, J. Kwiecinski, A.D. Martin and P. Sutton, Phys. 
Lett.B335 (1994) 220, and P. Sutton private communication. 

65. M. Kuhlen, Phys. Lett. B382 (1996) 441. 

66. A.H. Mueller, Nucl. Phys. B (Proc. Suppl.) 18C (1990) 125; J. Pys. 
G17 (1991) 1443 and J. Bartels, A. de Roeck and M. Loewe, Z. Phys. 
C54 (1992) 635. 

67. E. Mirkes, private communciation. 

68. HI Collab., S. Aid et al., PA03-049 contr. to ICHEP’96, Warsaw, July 
1996. 

69. J. Bartels et al., Phys. Lett. B384 (1996) 300. 

70. HI Collab., S. Aid et al., Phys. Lett. B358 (1995) 412. 

71. HI Collab., S. Aid et al., DESY-96-215, Hamburg, 1996. 


24 





ZEUS 1994 preliminary 

2 3 4 ln (Q) 



Figure 2: ln( l/x p ) max as a function of Q. The ZEUS (1993 & 1994) and Hf 
data are compared to results from OPAL, TASSO and TOPAZ. A straight 
line ht of the form ln( l/.r p ) maa . = bln(Q) + c to the ZEUS ln(l/.r p ) maa . values 
is indicated as well as the line corresponding to b = 1, discussed in the text. 
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Figure 1: Twice the measured charged particle multiplicity (for x ~ 10 -3 ) in 
the current region of the Breit frame as measured in DIS by this experiment 
and by the FNAL vp experiment at x ~ 0.1. The factor of two arises from 
the fact that there is a qq pair in e + e“ compared to a single quark in DIS. 
Also shown are the results from e + e“ experiments. 
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Figure 3: Preliminary inclusive scaled momentum cross section in the Breit 
current region for ep —> eX from this experiment. The different symbols 
represent the different .r-bins. 


10 




